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Abstract 

By considering a lattice model of extended phase space, and using techniques of noncommutative differential 
geometry, we are led to: (a) the conception of vector fields as generators of motion and transition proba- 
bility distributions on the lattice; (b) the emergence of the time direction on the basis of the encoding of 
probabilities in the lattice structure; (c) the general prescription for the observables' evolution in analogy 
with classical dynamics. We show that, in the limit of a continuous description, these results lead to the 
time evolution of observables in terms of (the adjoint of) generalized Fokker-Planck equations having: (1) a 
diffusion coefficient given by the limit of the correlation matrix of the lattice coordinates with respect to the 
probability distribution associated with the generator of motion; (2) a drift term given by the microscopic 
average of the dynamical equations in the present context. These results are applied to fD and 2D problems. 
Specifically, we derive: (1) The equations of diffusion, Smoluchowski and Fokker-Planck in velocity space, 
thus indicating the way random walk models are incorporated in the present context; (II) Kramers' equation, 
by further assuming that, motion is deterministic in coordinate space. 
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1 Introduction 



In the last few years there has been an increasing interest in exploring the possible relevance 
of noncommutative geometry (NCG) in various areas of physics. Roughly speaking, the 
basic idea in this context, is to try to develop an appropriate conceptual and mathematical 
framework in which the fundamental object is no longer a manifold (intuitively conceived as 
the generalization of a geometrical surface), but rather an appropriate algebra A (e.g. of R 
or C functions on it). This is (partly) motivated by the well-known result that commutative 
C*-algebras correspond biuniquely to locally compact topological spaces, hence that all rel- 
evant information for the topological structure of such spaces is encoded in their algebra of 
functions (e.g. ]]]). 

From an algebraic point of view there are two possible lines of approach, depending on 
the commutativity or noncommutativity of the algebras considered. Though the latter has 
received considerable attention, especially in connection with the study of quantum groups 
and quantum field theory (|2j, [|, f|, |j, |j and references therein), it is true that even if the 
algebra A is commutative, the resulting structures are quite rich, allowing for a considerable 
generalization of familiar geometric concepts, while keeping a rather close contact with ordi- 
nary (commutative) differential geometry. In particular, geometrical concepts and methods 
may be developed on both "continuous" and discrete sets, with interesting applications in 
such diverse fields as riemannian geometry, gauge field theory, integrable dynamical sys- 
tems, stochastic calculus, symplectic mechanics and kinetic theory (0-0]; for a review 
HJ and references there in). In this connection, a basic tool is provided by the 



sec 



(great variety of) differential calculi that can be defined on a commutative algebra A with 
unit 1, generalizing the familiar (exterior) calculus of differential forms. More specifically, a 
differential calculus (DC) on a commutative algebra A, is a graded associative, in general 
noncommutative algebra 

n(A) := n r (A) 

reNo 

where Q°(A) = A, Q r (A) is an ^4-bimodule, equipped with an exterior derivative, i.e. a 
linear operator 

d-.n-^n, d 2 = 0, d(l) = (1.1) 



d(uuj') = (doj)u' + (-l) r uJ da/, u G Vl r (A) (1.2) 

Elements of Q r (A) are called (differential) r-forms and we assume that as an ^4-bimodule, 
Q r (A) is generated by d^ 1 ^)). 

A key concept here, allowing for interesting conceptual insights, is the commutative and 
associative product of 1-forms defined byQ 

iW/:= [u,f] : = uf-fu, fEA,ue Q 1 (1.3) 

1 From now on we write fl r instead of fl r (A) and the symbol ":=" indicates a defining relation for the 
object that lies on the side of 
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and extended bilinearly in A, so that 

(fug) • [f'u'g') = ff\u • u')gg', f, f, g, g' E A, u, J e fi 1 (1.4) 

(for more details, see []2"0| , §§2, 3). For the existence of this • product, commutativity of A 
is essential. 

In the usual DC, eq( |1.3| ) is zero so that • is trivial. In fact, by ( |1.3| ), the Leibniz rule 
( |1 . 2|) can be written in Q l as 

d(fg) = fdg + gdf + df.dg (1.5) 

Thus (|1.5| ) stresses the nature of noncommutative DC as a deformation of the usual DC. 
On the other hand, it is formally identical with the generalized Ito differentiation rule in 
Stochastic Calculus (StC), /, g being semimartingales and the 3rd term in ( |1.5| ) is the so- 
called bracket of the processes, related to their quadratic variation (pl|, appendix eq(4)). In 
fact, it can be shown that this is not merely a formal analogy; more precisely, a noncommu- 
tative product can be easily defined for semimartingales and their differentials, so that the 
stochastic differentiation rule becomes identical to ( |1 . 5|) (|| §3). In this way one can obtain 
results of StC by employing techniques of NCG (e.g. see [22| §3 for the Ornstein-Uhlenbeck 
process) 

Motivated by the above remarks, one may develop differential geometry on the basis of 
a minimal deformation of the usual DC, namely, a calculus in which 

df»dg*dh=[[dh,g],f] = (1.6) 

and explore its relation to StC. Further motivation for this, stems from the following remarks: 
From a physical point of view, StC is a mathematical formalism on the basis of which 
a precise meaning can be given to stochastic dynamical models of physical systems whose 
time evolution cannot be given in terms of deterministic flows in their phase space (e.g. 
brownian motion). These are systems for which the concept of a trajectory in phase space 
is not defined. Mathematically speaking, their evolution cannot be described by semigroups 
of Perron-Frobenius operators, since their states often evolve under a 2nd order differential 
operator (differential operators can generate Perron-Frobenius semigroups only if they are of 



the first order, see e.g. [p3| , ch.7). The appearance of such operators is a basic feature of StC 
(via Ito's formula). In this way, by appropriate generalizations of basic geometrical concepts 
in the context of StC, a general mathematical framework results for the description of such 
nondeterministic systems ([0])- On the other hand, 2nd order evolution equations are at the 
heart of kinetic theory, for describing time-irreversible evolution towards equilibrium; they 
are obtained by employing some approximation scheme to the time-reversible microscopic 
dynamics. 

Taking account of the above discussion and starting with (|1.6|) , we have formulated basic 
concepts of differential geometry, tensor calculus and symplectic mechanics in a noncommu- 
tative context and we have shown that hamiltonian equations for observables are (adjoint 
to) generalized Fokker-Planck (FP) equations like those encountered in Ito's StC and the 
kinetic theory of open systems (@, fll5"l ). 

These are suggestive mathematical developments for the relevance of NCG to StC and 
kinetic theory. However, from a conceptual point of view, the physical meaning of noncom- 
mutativity has to be made clearer, particularly in connection with the time evolution of 
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physical systems. On this issue, a basic intuitive idea is that noncommutativity is due to the 
fact that differentials (i.e. 1-forms) have a "size", so that different results are obtained when 
functions are evaluated at the left end of the 1-form (its multiplication by functions from 
the left) and at the right end of it (multiplication from the right) - cf . Hl3|| , §6. On the other 



hand, 2nd order FP equations result in the context of kinetic theory (based either on micro- 
scopic dynamics or stochastic models), by assuming that, although individual microscopic 
interactions happen on a time-scale very short compared to the (macroscopic) relaxation 
time of the system as a whole, yet this microscopic time-scale is long enough to allow for 



many changes in the microscopic configurations of the system (see e.g. [24| §§II.l, II. 4 



ch.IX.4). Therefore, although eventually one passes to the limit of a continuous description 
in terms of macroscopic kinetic equations, the fine structure of the system on the "infinites- 
imal" microscopic time-scale has already been taken into account. This is a methodology 
widely used in statistical mechanics and kinetic theory; to start from a discretized picture of 
the system and subsequently pass to a continuous macroscopic description in some appro- 
priate limit (e.g. thermodynamic limit, appropriate scaling of the physical quantities, |26) 
ch.3.3, 0, §1, (28|, chill). 

On the basis of this, we may try to develop a discrete analogue of the formalism in 
|i"3fl , [15| , aiming at the following: (a) To throw more light into the physical meaning of 



noncommutativity of the DC; (b) to understand better the similarity of the DC to StC 
and (c) to provide a method for deriving irreversible evolution equations, which describe 
nondeterministic motions and which can be given a probabilistic interpretation, as exact 
rather than approximate results, much in the spirit of [I3j, §1, and in contrast with the 



general methodology underlying statistical mechanics (cf. the discussion in |13j §1).Q In 



this connection, the techniques developed in [O], [O for NCG on discrete sets will be used 



since, on the one hand they allow for a clear geometric representation of noncommutative 
DC and on the other hand, they are conceptually closer to those used in the study of specific 
models in stochastic dynamics or statistical mechanics.^ 

In view of the above, in the present paper we will show how simple model kinetic equa- 
tions can be derived in the context of DC on a discrete set, having a three-fold aim: 

(i) From the conceptual point of view, to clarify the relevance of noncommutative DC to 
kinetic theory and more generally, to the description of physical systems in terms of irre- 
versible evolution equations; 

(ii) since in kinetic theory and StC, probabilistic concepts are central, to explain how such 
concepts could naturally arise in the context of DC applied to systems conventionally studied 
by kinetic theory and stochastic dynamics; 

(iii) in analogy with classical dynamics, to formulate a general prescription for what dynam- 
ical evolution means in a formalism based on NCG. 

More specifically, the paper is organized as follows: 
The basic ideas in this section and the above general aims of the paper, are already appar- 

2 Notice that, implicit to the derivation of kinetic equations by approximating microscopic dynamics, is 
the assumption that the higher the approximation, the better are the resulting equations. However, this 
usually leads to differential equations of order higher than the second, in direct conflict with the fact that 
at least for linear autonomous equations, it is known that only differential operators of order at most 2 can 
generate semigroups of solutions admitting a probabilistic interpretation, fl30| . 

3 Though formalisms based on the "continuum" , more often than not can be manipulated mathematically 
more easily than formalisms based on discrete methods, the latter often provide deeper conceptual insights. 
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ent in the ID-model considered in section 2. This fact suggests the line of approach in the 
subsequent sections. The basic ideas are introduced systematically in section 3. Specifically, 

(a) the extended phase space as a discrete set, in particular a hypercubic oriented lattice, 

(b) vector fields as generators of motion and probability distributions, (c) the emergence of 
the time direction on the basis of the encoding of probabilities in the lattice structure of the 
extended phase space, (d) the general prescription of the evolution equation for observables, 
in direct analogy with conventional dynamics. In section 4 we apply the previous ideas and 
results in ID-problems, thus getting the incorporation of random walk models in our context 
and accordingly showing how the diffusion, Smoluchowski and FP equation in velocity space, 
result in the continuous limit. In section 5 the approach of section 4 is generalized in N di- 
mensions, giving in the continuous limit a generalized FP equation in which, (1) the diffusion 
matrix is nonnegative-definite and is the limit of the correlation matrix of the lattice, with 
respect to the probability distribution associated with the vector field which is the generator 
of motion; (2) the drift terms are specified by the dynamics of the particular problem under 
investigation. Some comments on the form and properties of the transformations in phase 
space, are given in section 6. In section 7 the results of section 5 are applied to the 2D-case. 
As an example, Kramers' equation is derived by assuming Newtonian equations of motion 
with friction, and requiring that in the continuous limit, motion in coordinate space is deter- 
ministic (i.e. trajectories exist there). In section 8 the nature and mathematical properties 
of the limiting procedure used in this paper and of other such procedures, are examined, 
showing that only the limit used here is well defined for all possible (1st order) differential 
structures on the lattice. Moreover, if the conceptual framework of §3 is employed, then we 
show that the only differential operators which can be generators of evolution of observables 
are at most of order 2. Finally, in section 9 we summarize our results and discuss the main 
ideas involved, as well as possible further elaborations in the present context. 

2 A One-Dimensional Model: Motivation For Further 
Developments 

As already mentioned in §1, there is a formal similarity of the Leibniz rule in noncommutative 
DC, eq (|1.5| ), and stochastic differentiation. Actually, in |§, it was shown that Ito's stochastic 
differentiation can be made an exterior derivative in noncommutative DC. In this section we 
explore further this issue by means of a simple ID discrete model, following the rationale 
of §1. Specifically, we consider a point moving in one dimension x, so that motion can 
be described in M 2 with local coordinates (t, x), equipped with a DC having the following 
commutation relations (t is the time; for the choice of the signs see §3.4) 

dt»dt = -bdt, dx»dx = -hdt, dt»dx = -bdx, h,b>0 (2.1) 

and we assume that {dt, dx} form a basis of 1-forms. Associativity of • implies that 

dx • dx • dx = a 2 dx, a 2 := hb (2.2) 

so that in the limit a, b — > 0, a? jb — > h = constant, we recover in one dimension, the DC of 
H §5, |]13| §4, in close analogy with the Ito StC and with commutation relations 

dt • dt = dt • dx = 0, dx»dx = -h dt (2.3) 
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It is a special case of the DC defined by With the transformation 

1 fx t\ 1 fx t\ . , , . . . , 

u = -\ ) , v = -- - + r & t = -b(u + v), x = a(u - v) (2.4) 



2 \a bj 2\a b / 

( p.l| ) becomes 

du • dw = dw, dw • df = 0, dv • dv = dv (2.5) 
which is the DC on an oriented square lattice. Hence, from ( |1.3| ) we easily get 

d/ = (f(u + l,v)- f(u, v)) du + (f(u, v + 1)- f(u, v)) dv (2.6) 
so that a direct calculation using ( |2.6| ) gives 

df{t, x) = B x f(t - b, x)dx + (d- t f(t, x) - ^A x f{t - b, x)) dt (2.7) 



where 



(8Mt,,):= f{t ' X + a) ;J {t ' X - a) , (2.8) 



;= /(«■») -/(«->■*) (2.9) 



(A x /)(t,x) := /(^ + °) + /(^ - °) - (210) 

a 2 



In the limit a, 6 — >• 0, a jb — > h = constant, ( |2.7| ) becomes 

d/(*,s) = 9 s /(t,a;)dx+ ^/(t,x) - ^/(t,x) ) d/ (2.1 I ) 



recovering in ID the results of previous works (|T3[ eq(4.13), eq(5.3), [^] eq(5.1)) and 
essentially Ito's formula for the differential of a function of a Wiener process with variance 
h. 

In order to explore this similarity further, we consider the (free) motion of a point in this 
context, by introducing the concept of a differentiable motion in analogy with the ordinary 
DC. To this end, let (M, s, d') be a discrete DC on R (s being a coordinate function on M 
parameterizing this motion), defined by 

d's.d's = -bd's (2.12) 

Then, we define a differentiable motion to be a mapping 

7:M^M 2 , 7 (s) = (t(s),x{s)) (2.13) 
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inducing a "pull-back" homomorphism 7* on fi 1 (IR 2 ,d) to d') withQ 

%(x)=x(s), 7,(t)=t(s), 7, od = d'o 7 , (2.14) 

Eq( j2.13| ) implies 

/ , \ x(s) — x(s — b) u . , . . , . t(s) — t(s — b) ,. ... 
7,(dx) = — ^ - d's =: z(s)d'a, 7*(dt) = — — ^ d's =: t(s)d's (2.15) 

Applying 7* to the commutation relations ( |2.1| ), ( j2.2| ), we get the following consistency 
relations 

bx 2 {s) = hi{s), bx{s)i{s) =bx{s), bi 2 (s) = bi{s) (2.16) 
Therefore, with x(s) 7^ 0, £(s) 7^ to avoid trivial cases, we get t = 1, x = ±a/b, hence 

x(s) -x(s-b) _ a 

b ~ b [Z - U) 

The following remarks can be made here: 

(i) Motion consists of jumps ±a in discrete time intervals b, i.e fl2.17Q gives the discrete 
velocity of the point. In the continuous limit which leads to Ito's StC, (|2.17|) becomes 
infinite, strongly suggesting the picture of brownian motion as a random walk with equal 
left and right probabilities of jumps ±a at time intervals b and diffusion constant a 2 /b (see 
e.g. @ eq(17), 0, eq(7)). 

(ii) The DC defined by Q2.1|) is up to a coordinate transformation the DC on an oriented 
square lattice (cf. fl2.4Q ). 

If the model is to be taken seriously, then these remarks raise the following issues: 

• How can probabilistic concepts be introduced, given that no such concepts appeared 
above? In other words, is there some deeper reason for the formal similarity of the 
model with random walk, which may lead to the introduction of probabilities at a 
fundamental level? 

• The above model corresponds to ID-free motion. How can it be generalized to higher 
dimensions and/or other models'? That is, what is the general prescription of time 
evolution in the context of noncommutative DC? 

These questions correspond to the general aims of the paper, expressed in §1 and are discussed 
in the next section. 



3 Basic Ideas And General Formalism 

In §1 we noticed that in statistical mechanics and kinetic theory, one often starts from a 
discrete picture of the system under consideration, either for methodological reasons, or 

4 More generally, a mapping <f> : (Q'(A),d') — > (f2(.4),d) between two DC, is called differentiable, if it is a 
homomorphism of Q'(A) to Q(A) and cf> o d' — d o cf> (see |12], §4). 
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because at the microscopic level these systems are in some fundamental sense discrete (e.g. 
composed by a number of particles). This idea also appears in other areas, like lattice 



field theory (e.g. f32fl , ch.8), or certain approaches to quantum gravity (e.g. |33|, P4[ , see 
also ||35|| , [|36f| for recent reviews, and references in [0 §1). Therefore, motivated by the 
discussion in §§1 and 2, we start with a discrete (finite or denumerable) set A4, which we 
call the "phase space" and eventually pass to the special case of an oriented hypercubic 
lattice Z N+1 . We will use the concepts and techniques developed in |Tl|, [12], [L7|, a summary 
of which is given below. At this point we stress that, as explained in §3.1, starting with 
a discrete picture has far-reaching consequences, since any differential structure on A4 is 
necessarily noncommutative, i.e. discretization implies noncommutativity of the DC on M.. 
Thus, the need for, or necessity of, a discrete description of physical objects, can be seen as 
a basic motivation (or explanation!) for introducing noncommutative differential structures 
in the study of physical systems. 

3.1 Differential calculi on a discrete set: A summary 

Let A be the algebra of C- valued functions on A4, with the usual algebraic operations. Then, 

f=J2f(i)ei=:J2f^ W f eA M./): 'V ^ G A (3.1) 



i 

We introduce a DC on A as in §1. By defining 

-••= {•■;•• Hi (3-3) 

so that 

f &ij ft f fj &ij (3-4) 

we get ejdej = — Yuj^i e ij an d consequently 

de i = J2( e H ~ e ij) ( 3 - 5 ) 

3 

On the basis of this we can prove that {e^} is linearly independent and a basis of Q 1 taken 
as a vector space over C. In particular 

d/ = y~](fj - fi)eij , df • e i:j = (fj - fye^ , ey • e H = 6^6^ , V/ G A (3.6) 

h3 

It can be proved (|lT| §11) that the induce a basis on Q r , r > 2 over C via 

"{cii---'i r . I ^-ii—ir ■ €-i 1 i 2 ei 2 i 3 ■ ■ • — j^,.} (3-7) 
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Therefore, there is a simple way to obtain any particular DC on A simply by imposing 
relations among the ey, which turn to be equivalent to putting some of the equal to 0. 
This implies additional form relations for r > 2. Actually, the DC defined by ( |3.3| ), ( |3.7| ) 
is the largest one, called the universal DC on A. At this point it is conceptually suggestive 
to notice that with i ^ j may be represented as an arrow from i to j, i >— > j. Then, 
universality means that all pairs of points of M. are connected by two antiparralel arrows 
and any DC on A is obtained by simply discarding some arrows from this complete di-graph 
(viz. directed graph). In this representation, ( |3.5|) gives 

dej = (sum of incoming arrows at €)— (sum of outgoing arrows at i) 

and similar interpretations hold for ( |3.7| ). Moreover, fl3.4| ) shows that left (right) multipli- 
cation of 1-forms ejj by functions, implies that functions are evaluated at the starting (end) 
point of the corresponding arrow. This proves that any nontrivial DC on a discrete set M. is 
necessarily noncommutative, since at least one is nonzero so that d/«ejj is not identically 
zero. Actually, this quantity is (proportional to) the change of / along the arrow i >— > j. 
These remarks are important for what follows in this section. 



3.2 Vector fields on a discrete set M. 

Vector fields X on A are elements of the ^4-bimodule X(A) dual to Q , defined by 

(feu, gX ■ h) := fg(uh , X) , /, g, h E A, u E Q 1 (3i 

where ( , ) denotes duality contraction with Q 1 considered as a left ^4-module. 
By (|3.8j), vector fields act as operators on A 



X(f):=(df,X) (3.9) 

hence by (|3.8j ) and the Leibniz rule 

(fX)(g) = fX(g) , (X • f)(g) = X(fg) - gX(f) (3.10) 

We define elements {dy} G X(A), "dual" to {e^}, by 

(e%j,dki) = Sikdjiet (3.11) 
(putting dkk = 0). A direct calculation, using (|3.11|) , (|3.4j) , gives 

ekdij = Sikdij , dij ■ et = Sjkdij (3-12) 

so that 



<h r f fUVhj- f^i !(<)»,,■ d ij {f) = {f{j)-f{i))e i (3.13) 
Hence for any X G X 

X = ]TX%, X{f) = Y,ih-h)X ij e i , XV EC (3.14) 



i,3 i,3 
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and we may write X as an ordinary difference operator X = ^\ . X^dij =: Y2 a X a d a , where 
a is an index for the arrows of the di-graph. The following proposition plays an important 
role in this paper and its proof is straightforward: 

Proposition 3.1 Let I be the identity mapping on A. Then <fi : A — > A is an endomorphism 
of A, if and only if 

X:=(f)-I (3.15) 

satisfies 

X{fg) = gX{f) +fX{g) + X{f) X{g) (3.16) 

We remark here that X in general is not a vector field on Ai . Nevertheless eq( |3.13|) readily 
implies that 9y (equivalently, d a ), satisfy ( |3.16|) . Using this, we easily prove: 



Proposition 3.2 A vector field X G X(A) generates an endomorphism <p via ( \3.1b] ), if and 
only if 

X = J2 x<xd * , with X a X 13 - 8 af} X a = (3.17) 

a 

Eq( |3.17|) simply means that at each point of Ai, at most one of the X a is nonzero and equal 
to 1. By ( |3.16j ), ( p,17| ), X is the generator of a discrete semigroup of endomorphisms of A, 
{<j) n = (I + X) n , n G N}. 

As it is shown in appendix A. 1, is an automorphism, if and only if it induces a 1-1 and 
onto mapping $: Ai — > Ai. In this case, is essentially the Koopman operator associated 
with $ and maps the basis {e;| i G Ai} onto itself (see e.g. [^3| §3.3). Therefore, the orbits 
of the semigroup generated by X in ( p. 17] ) are determined by the trajectories of points of 
Ai under $ (i.e. the flow defined by $), which may be thought as lying along arrows of the 
di-graph corresponding to the DC on Ai (notice that by the remark above, at each point 
of Ai at most one term in ( |3.17| i) is nonzero, and by ( |3.13 ) each dij gives the change of 



functions along the arrow i >— > j). Thus, vector fields generating automorphisms of A, define 
flows on Ai along arrows of the di-graph that map points of Ai to points of Ai. In this 
picture, stemming from the discrete character of Ai, motion along a trajectory defined by 
a vector field X, is constrained to be along arrows of the di-graph, so that at each point i 
of Ai (i.e. vertex of the di-graph), motion along any arrow is either certain or impossible, 
depending on X. 

The above discussion raises naturally the question what kind of "motion" is described by 
an arbitrary vector field, i.e. one not necessarily satisfying 17[ )? Evidently, {e^} (equiv- 
alently, Ai) is no longer mapped onto itself, but to a set of nontrivial linear combinations 
of the e^'s, that is, to a "superposition of points" (cf.0, §3). This reminds the situation 
in statistical mechanics; non-unitary semigroups of operators on the state space, transform 
pure states into mixtures in the quantum case and ^-distributions into more general ones in 
the classical case. 

On the basis of this analogy and the discussion in the previous paragraph, we shall inter- 
pret a vector field X as the generator of evolution on Ai, having at each point a particular 
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well-defined probability associated with each arrow that emerges from that point. More pre- 
cisely, in this interpretation, the component X l i of X gives the transition probability for the 
"infinitesimal" motion from i to j. 

In this picture, vector fields X acquire a double role; as generators of evolution and 
as states giving the transition probabilities for "infinitesimal" changes. This double role is 
reflected in the notation of ( |3.9| ) so that (df, X) = X(f) is the expectation (average) value for 
the "infinitesimal" change of the observable / along X (cf. fl3.14|) ). It is a key idea following 
from adopting a discrete structure as our starting point. More precisely, though both in 
the discrete and the usual continuous case, vector fields are linear combinations of 1st order 
difference and differential operators, respectively, it is only in the continuous case that any 
such combination generates automorphisms of the algebra of functions, or equivalently, flows 
of trajectories (i.e. what we call here a deterministic evolution, or motion). This is due to 
the fact that, in contrast to the continuum, nontrivial linear combinations of displacements 
do not give a permissible displacement on the di-graph. In the rest of this section we shall 
pursue further this idea. 

3.3 The oriented hypercubic lattice and vector fields on it 

To be more specific and for applications in subsequent sections, we shall henceforth take M. 
to be an oriented hypercubic lattice; that is, M. = 7L N with elements k e M. and 






{/t} being the canonical basis of K. 



(for details see |TT|). Then, by defining 




(3.19) 



k 



and writing u := (u M ), we get for f E A 




(3.20) 



k 



so that from ( |3.6|) 



(13T4D we get 





k 



df 





'+-UM — 



k ^k,k+fi 



(3.23) 
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where we use the same symbol for / as a function of k and of u, to avoid a cumbersome 
notation. 

Clearly, {du^} is a basis for 1-forms with dual basis for vector fields {<9 +m m} 

x = £ ^ =: J2 P"d + u, , F" = J2 X%U ^ ( 3 - 24 ) 

jfc,jj A* jfe 

(du", = 6£, du" • du" = d^du" (3.25) 

the last relation following from ( |3.6|) . 

The interpretation of vector fields as states giving transition probabilities, discussed in 
the previous subsection is manifested in ( 3.24 ); X k,k+ ^ is the transition probability at a point 



k in the direction of the axis u M , so that P M is the distribution of this probability on Ai. 
Therefore, from now on we consider vector fields X such that 

P^>0, ^P" = l (3.26) 

Thus, in view of (|3.22|) , (|3.25|) , (df, X) gives at each point of A4, the expectation of df (the 
"infinitesimal" change of f), along X , with respect to the probability distribution determined 
by X. To see that this interpretation is consistent and for later use as well, we notice the 
following: The l.h.s. of ( |3.17| ii) is rewritten for an arbitrary X = P m <9 +u m as 



P ^ ;= S r>pn _ pupv = ( dM M . du v } X ) - (du", X) (du u , X), (3.27) 

where ( p.25| ) has been used. (P^) is the correlation matrix of the ^-dimensional random 
variable (dw M ) having a probability distribution (determined by) X, where • is the natural 
(intrinsic) product of 1-forms. P^ v vanishes if and only if X generates a flow of trajectories in 
M.. Actually, at k e A4, P^ u {k) is the correlation matrix of the following random variables 

j (t) — / 1 ^ change at k is along the u M axis 
M \ otherwise 

with probabilities P"(jfc) (cf. 0, §3.3, p.331). 

The following proposition has an interesting probabilistic interpretation to be used later 

on 

Proposition 3.3 For uj = ^^s^du^ e Q 1 , write s t = (s M ) (row vector). Then 

(umu,X) - ((^,X)) 2 = s'Ps (3.28) 
with P = (P M "). Moreover, for 

p:=J2 dufl = J2%k + » ( 3 - 29 ) 
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we have 



p • uj 



(3.30) 



(pm Pi X) - {{p,X)f = 0, (p»uj,X) - (p,X)(uj,X) = 0, WujEQ 1 (3.31) 



= ((du" - (du", A)p) • (dw^ - (du u , X)p) , X) (3.32) 

The proof is immediate, but we stress the fact that, for (|3.31|) it is essential that P M = 1. 
Eq (|3.32| ) is an equivalent expression of (|3.27| ), obtained by using the unit p with respect 



to the natural multiplication • of 1-forms. On the other hand, in the present conceptual 
framework, ( |3.28| ) gives the variance of the random variable (s M ) (and for that matter, of uj) 
having probability distribution X. Thus ([3.31 ) says that the unit of (Q, •) has zero variance 



and is uncorrelated with all uj e Q 1 . Thus, as a random variable, p is constant, a fact that 
will be used in the next subsection. Algebraically, ( |3.31| ii) says that P has always a zero 
eigenvalue and a corresponding lD-eigenspace {A(l, 1, . . . , 1) , A G R}. It is easily seen that 
this is the only generic eigenspace of zero, in the sense that it is the only one existing at all 
points of M. 



3.4 The concept of time 

As already mentioned at the beginning of this section, we consider M. to be the extended 
phase space, in the sense that it contains the "time axis" . Adopting a newtonian picture we 
require: 

(a) For every evolution (e.g. motion) in Ai, a change of time is required; 

(b) Time flows with certainty, that is, there is always a change in the time axis. 

Then, if any evolution on M. is determined by some vector field X, then (a) means that the 
time change is given by 1-forms r such that (r,d +u ^) ^ at all points of M.. Moreover, in 
view of the discussion following proposition 3.3, (b) implies that 

T = bp = d(-&J]V) =: dt, (dt,X) = -b, dt»dt = -bdt, b 6 R + (3.33) 
so that (a) is also satisfied. 

(i) Since VtF" , P^ are dimensionless, b is put for dimensional reasons. 

(ii) By the interpretation of vector fields X both as generators of evolution and states, 
X(f) = should describe the evolution of observables in the extended phase space. Then, 
as explained in appendix A. 2, X(t) should be negative, in analogy with classical dynamics, 
hence the choice of the (— ) sign in ( p. 33 ). 



3.5 General prescription of dynamical evolution 

In summary, the approach described in this section leads to a general prescription for time 
evolution in j\A: 
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M. is taken to be the extended phase space (cf. §3.4). 



• Via ( |3.9| ), vector fields X are seen both as generators of evolution and as states giving 
the transition probabilities for "infinitesimal" changes (cf. §§3.2, 3.3). 

• Time change appears in steps of duration b. 

Therefore, elements / of A are observables, whose "infinitesimal" change has an expectation 
value with respect to X, (df, X) = X(f). But then, by remark (ii) in §3.4 and appendix A. 2, 
—X(f)/b gives the rate of change of f along X in extended phase space, hence, in analogy 
with classical dynamics, we take as the general dynamical evolution equation 

X(f) 

evolution equation for observable /: — = (3.34) 



which becomes with the aid of (|3.22p , (|3.24j ) 

- 1 E p "( u )(f( u + £) - /(«)) = ( 3 - 35 ) 

Eq( |3.34j ) is a general prescription, independent of the choice of M. as an oriented hypercubic 
lattice. Nevertheless, it is our starting point in the applications in the next sections. As we 
shall see there, in the continuous limit, ( [3.35D reduces in particular cases to the adjoint of 
well known kinetic equations, i.e. to evolution equations for observables. 



4 One Dimensional Problems 

It is quite straightforward to apply the general ideas and results of §3 to specific cases. 
To keep technical details to a minimum, we consider in this section as an illustration, the 
derivation of well known ID-model evolution equations and leave the general treatment to 
the next sections. 

Motivated by §3.4 and the model of §2, we consider a coordinate system (t,x) in the 
notation of §2. Although the transformation from the lattice axes to (t, x) can be any 
invertible transformation, we postpone this till the next section and use "light-cone" coor- 
dinates defined by ( |2.4j ). Then, for / e A, df is given by (|2~7|). On the other hand, ( |3.24j ) 
is 

X = pd +u + qd +v = (du, X)d +U + (df , X)d +V 

y k J \ k ) 

where p, q > and p + q = 1. Using fl2.7Q in ( |3.34| ) and that p + q = 1, we readily obtain 
the evolution equation 

1 d- t f(t, x) - -AJ(t -b,x) + d x f{t - b, x)-^ 1 

d-tfit, x) - ^AJ(t - b, x) \ - ~(p - q)d x f(t -b,x)=0 (4.2) 
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In the derivation of ( |4.2| ), we remark that the following features are present in any number 
of dimensions (see e.g. §5 and appendix B): 

(a) Because of flOf]), (dt,X) = -b. 

(b) 2nd order difference operators appear only in the coefficient of dt. 

(c) the term involving (dx, X) = X(x) has the form of a drift term. 

More precisely, (dx, X) has been interpreted in §3.5 as the "infinitesimal" change of x in one 
step of time b, hence it is natural to assume that it is proportional to b. This assumption can 
be justified on the basis of the general approach described in §5 (see ( |5.13| ), the derivation 
of ( |5.15| ) and appendix A. 3). Since we eventually pass to the continuous limit considered in 
§2, namely a, b — > 0, a 2 /b — > h — constant, we have 



and therefore 



(dx , X) = a(p — q) = o(b) =>- p — q = o( — ) = o(a) (4.3) 



(dx,X) a, . , A tS 

\_I_/=_(p_ g )_>i2 (4.4) 

say, when a, b — > 0, a 2 /b — > h, so that in this limit, ( |4.2j ) becomes 

d t f - RdJ - h 2 J = (4.5) 

To see the physical meaning of this evolution equation, we go back to ( |3.24|) and consider 
two special cases. 

1. x^ Ui ,X^ %+2 are constant: 
Then from Q, (Q) 

x U+i = p = . 1 _ ia ; x fc,fc+2 = 5 = . 1 + 7a ? 7 = constant (4. 6 ) 

As it is evident from the figure, p, q represent the transition probability for right and left 
motion respectively, along the x-axis. This is a random walk model in ID-position space 
and in view of (fO|), the evolution equation in the continuous limit, is 

dtf + 2jhd x f - h 2 J = (4.7) 

This is the adjoint of Smoluchowski's equation for a constant force field, of intensity propor- 
tional to 7 and diffusion constant h, derived here in a conceptual framework totally different 
from that conventionally used (e.g. the theory of Markov processes; see [^T[, §3). 

2. Evolution in velocity space. 

Here x represents the velocity of a particle and we take (cf. ( |4.3j )) 

X % ' Ui = -- (3b(k l - k 2 ), X^ +2 = - + pbik 1 - k 2 ) , (3 = constant > (4.8) 

hence by ( ggg ), (fZ|) 

11 b 

p= (3b(u -v), q = - + fibiu - v) , p- q = -2-f3x (4.9) 

2 2 a 
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Coordinate transformation on the square lattice 



so that in the continuous limit, ( |4.2|) gives 

d t f + 2(3xd x f-^dlf = (4.10) 

This is the adjoint to the ID FP equation in velocity space, giving the Ornstein-Uhlenbeck 
process (diffusion constant h, drift coefficient (3). We may remark however, that if x is 
interpreted as a position coordinate, then ( |4.10| ) is the adjoint to the Smoluchowski equation 
for a harmonic external field (cf. ||31|| , §4). 

5 Evolution Equation: The General Case 

In this section we apply the results of §3 in the case of an (N + l)-oriented hypercubic lattice 
M. 

The lattice coordinates are dimensionless. We consider linear coordinate transforma- 
tions to physical coordinates 

/ = «^^, n,u = 0,l,. ..,N (5.1) 

V 

in which are scaling parameters with the dimensions of x M . In view of §3.4 we require 
that 

x° = t, hence a = -b, A° = 1 (5.2) 

and we shall eventually consider the continuous limit 

a„-0, f^ ±00 ' (5.3) 

with the depending on the a^'s, in general. Eq( |5.3| ) is motivated by the limiting procedure 
considered in ]13|, [T5| and §2 here. Nevertheless, its nature is further explored in §8. Here 
we only notice that, in contrast to (t,x l ), the w M are not defined in this limit. Moreover, to 
avoid a too technical development, we restrict the discussion to linear transformations with 
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constant coefficients, thus making the physical and conceptual aspects of our approach more 
clear. In what follows, latin indices run from 1 to N, greek ones from to iV and symbols 
with a hat denote quantities in the limit (|57 



By (5J.) we can write 

Qj v 

V V 

so that dH), (|3^5| ) become with the aid of Q 

df = (d- t f(t, x) - ±Af(t - b, x)j dt + J2 9if(t ~ b, x) dx l 



(5.4) 



where 



dif(t,x) 
Af(t,x) 



(5.5) 

(5.6) 
(5.7) 



and 



dt • dt = -b dt, dt • dx i = -b dx i 

dx' • dx j = -^p- Y KMB£ dt + Y ^±A\A{B^ dx r < 



(5. 



(1,771 



With the interpretation of vector fields X given in §§3.2, 3.3, we consider those X that 
satisfy (|3.26|) , i.e. 



X = ^P^ +u m, P^>0, J2 P>1 = 1 ( 5 - 9 ) 

hence (dt,X) = —6, (du\X) = P l . Therefore the evolution equation ( |3.34|) becomes (c.f. 



\x(f) = d_ t f(t, x) - l -Af(t - b, x) - £ T( AP y W -b,x) = (5-10) 



where, by Q5.1|) 



(dx\X) = at (APy 



(5.11; 



With the interpretation of vector fields introduced in §3, this is the average value of the 
"infinitesimal" change of x % along X in one time-step b. By assuming that / is sufficiently 
well-behaved function, Taylor expansion of (|5.10| ) around x 1 , gives 



\x(j) = - e t (apt ^ - ~ £ ^ fe w + 

i i,j \ /j, J 



o{a) = (5.12) 
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We proceed further in analogy with conventional dynamics. Let (R l ) be the generator of 
motion in this case (e.g. in Newton's or Hamilton's equations). Then, as it is explained in 
appendix A. 2, to the system of characteristic equations of the observables' evolution equation 
( |A.3| iii), are associated the 1-forms 

a 1 = dx i + ffdt (5.13) 

which vanish along the integral curves of the vector field that gives the dynamical evolution. 

Therefore, in the present discrete context, we postulate the existence of dynamics in the 
above sense, that is, the generator of dynamical evolution is determined by a similar relation, 
namely 

{a\X) = (dx\X) + R i (dt,X) = hence a^APf = bR { (5.14) 
(cf. (IAT5D . Hence, 

(APy = -R i = ^-R! + o(a) with R* -> (5.15) 

0{ tin 

Notice that the postulated existence of dynamics, that is (|5.14j i), implies that (dx l ,X) is 
o(b), thus justifying this assumption which was introduced in the derivation of ( |4.4| ). It is 
also equivalent to a similar relation determining the transition probabilities P M , as explained 
in appendix A. 3. 

Since (AP)° = 1 (because of (|5.2|) , (|0|)), we can invert ( |5.15| ) to get 

P" = Bj + ^^LB^r, hence in the limit = B% (5.16) 

rn 

so that by ( ^9|) 

= 1, BS>0 (5.17) 



Therefore ( |5.12j ) can be rewritten as 

- \x{f) = dt f-Y,^r h R i dif-\Y, V ij 9Af + o{a) = (5.18) 



Kb lJ 2 

hj 



where 

i7 y := ^£44 s o - r = h l3 ^2KK p " ( 5 - 19 ) 

In the continuous limit ( |5.3| ), we have 

dj - J2 R%f -\J2 ^ d Af =■ xf = (5-20) 
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This is our general evolution equation for observables /, corresponding to the DC obtained 
from (|5.8|) in the limit (|5.3|) 



dt • dt = 0, 



dt • dx l = 0, 



dx 1 • dx 3 



-rf dt 



(5.21] 



Eq( |5.21| ) is identical with the DC considered in |13|, [15| , called 2nd order DC there, because 
it leads to 2nd order evolution equations, i.e. ( |5.20 ). We shall come to this point again in 
§9. 

A number of interesting remarks can be made here: 

1) If ( J5.17D holds, then f\ %3 is a non-negative definite matrix, so that ( |5.20| ) is the adjoint of 
a generalized FP-equation with diffusion matrix rf-'and drift vector R % . 

2) Under the transformation (|5.1| ), the correlation matrix of the lattice, ( |3.27 ), becomes 



H 



(AFA 



(5.22) 



J2 a i a 3 A l A i P>1 - aiajiAPYiAPy = (dx i • dx 3 ,X) - (dx\ X) (dx 3 \ X) 

= (d(xV) - x l dx j - x 3 dx\ X) - (dx\ X) (dx j , X) (5.23) 
where we have used the Leibniz rule (|1.6|) . 



jjno _ q reflects the fact that the time change dt is the (essentially) unique eigenvector 
of zero of P, as discussed in 53.4 



Because of ( |5.15|) , the last term on the r.h.s of (|5.23|) is o(ctj 4 ). 
Using the unit p of the algebra (fl 1 , •), eq (|3.29|) , we may write 

H ij = ((dx* - (dx\X)p) • ((dx j - (dx j ,X)p) ,X) 



(cf. ( |3.32| )). In fact, H % i is the correlation matrix of the 1-forms a 1 , eq( ^.13| ), associated 
with the equations of motion via (|5.14|) 



H 



((a 1 - (a\X)p) . ((a j - (a\X)p) ,X) 
(^•a j ,X) - (a\X)(a j ,X) 
(a l »a j ,X) 



(5.24) 



Therefore in the continuous limit we get 



ffi = (dx 1 • dx 3 , X) = (d(xV) - x l dx J - x J dx l , X) 
= X(xV) - x i X(a: i ) - x^Xix^ 



(5.25) 



Thus, ff 3 is essentially the limiting value of (the matrix of) correlations of "infinitesimal" 
changes along the axes of the lattice, induced by the generator of evolution X. In fact, by 



( f5.22| ), ff 3 vanishes if X is the generator of a deterministic flow on the lattice, i.e. when the 
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concept of a trajectory is well-defined for X, as explained in §3.3 (cf. ( |3.17| )). This is reflected 
in the limiting equation ( |5.25| ), which gives the deviation of X from an ordinary derivation, 
that is, an ordinary vector field for which the concept of a trajectory is (locally) always 
meaningful. We may also compare rfi with the usual interpretation of the diffusion matrix 
in the kinetic theory of open systems. There, it is the correlation matrix of the hamiltonian 
vector field of the interaction hamiltonian of the open system with the bath, with respect 
to the state of the bath (see e.g. |38|] eq(4.16)). A similar interpretation of this quantity 
exists in the theory of markovian stochastic processes and of Ito's stochastic calculus (e.g. 
H chs.3, 4, [0, ch.VIII). 

On the other hand, the above comments provide an interpretation of the noncommuta- 
tivity of the DC defined by ( p.21| ) and studied in previous works (|13, TIJ): Starting with 



a discrete structure of the phase-space, somehow reflects the fact that originally one con- 
siders the fine microscopic structure of the system. At this level, changes are described by 
structures that are necessarily noncommutative and it is a generic characteristic of vector 
fields in this case, not to describe deterministic motion (i.e. flows of trajectories), in contrast 
with vector fields in the (ordinary) continuous case for which this is always the case. Then, 
the passage to the continuous limit corresponds to a coarse-graining process, in which the 
fine details of the micro-structure are no longer explicit, but which, nevertheless, have been 
taken into account. In particular, the above generic characteristic of vector fields on a dis- 
crete phase-space is expressed quantitatively in the continuous limit, by the non- vanishing 
of r) ij , eqflCTD- 



3) In 1) above, it was mentioned that (fj t3 ) is nonnegative-definite (being the limit of such 
matrices). On the other hand it is an easy matter to see that for such a matrix, fj u = for 
some fixed i implies ifi = for all j|] In our case, this means that with respect to i ! , X in 
( |5.20| ) is a derivation, hence it describes deterministic motion along the x*-axis. 

4) The limiting eq( |5.2U| ) is the most general linear autonomous differential (evolution) equa- 
tion, whose adjoint admits a probabilistic interpretation; more precisely, under mild reg- 
ularity conditions on the coefficients, it defines a positivity and normalization preserving, 
strongly continuous semigroup, globally on the space of (continuous) observables / having 
a finite limit at infinity. Its adjoint also defines a positivity and normalization preserving 
semigroup, satisfying a general if -theorem theorems 5.1, 6.1). This is assured by the 
non- negative definiteness of ifi . In our case this follows from ( |5.17|) , which is a consequence 
of the interpretation of a vector field as a transition probability-distribution for infinitesimal 
changes. Thus, if {d t , d x %} is the basis of vector fields dual to {dt, dx 1 }, then 

B% = (d< -b~d t ) = -bd t (u") - B% > 

Therefore, Bq > means that in the limit of the continuum (i.e. in a macroscopic descrip- 
tion — see 2) above), changes in phase-space are in the future direction (cf. the discussion 
following eq( |A.3| ) in Appendix A). Thus, in this picture, the " macroscopic" direction of 
time is intimately related to a probabilistic interpretation of evolution at the "microscopic" 
level (here the terms "macroscopic" and "microscopic" correspond to the terms "continuous 
description" and "discrete description" in the sense of 2) above). 

^{fj 1 ^) defines a (possibly degenerate) nonnegative-definite scalar product on M. N , and we apply the 
Schwartz inequality to the i-th element of the canonical basis of R with each j-th element j ^ i. 
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5) The expansion of the evolution equation fl3.34| ) (or ( 5 . 1 (J| ) ) in powers of the is not 



essential, but has been done because of the generality of the transformation (|5.1|). In each 



particular case X(f) may be written explicitly in terms of appropriate difference operators, 
which in the continuous limit yield the corresponding FP equation (e.g. as in §4). Such an 
iV-dimensional example is worked out in Appendix B. 



6 Transformations in phase-space 



There is considerable freedom in the choice of the coordinates (x l ), eq ([5.1|) , in such a way 
that the continuous limit exists. In this section we comment briefly on this issue. 

Let (x^), (x /fl ) be coordinate systems on A4, obtained from the lattice coordinates by a 
linear transformation like eq( |5.1| ) and let x' M = A^x u . In view of the discussion in §§3.3, 3.4, 
5, from a physical point of view we require that A leaves dt invariant. Then it has the form 



A 








A* 


A* 



(6.1) 



Such transformations form a group. Moreover, since dt is in the kernel of the correlation 
matrix (H tJ,u ), eq (|5.22|) , A leaves the form of (i.e. = 0) unaltered: 



H 















If we write (|5.1| ) in the form 



H' 











(AHA 



(6.2) 



[u 



A 





i— 1 


A 1 


4 



(6.3) 



where 1 = (1, 1, ... , 1), then 



H' = (AA)F(AAf 



(6.4) 



A natural (partial) fixing of the coordinate freedom is to choose A so that the correlation 
matrix, hence rfi in ( p.20| ), is diagonal. However many other possibilities exist, see e.g. 
Appendix B and section 7. 



7 Further Applications: 2D Evolution Equations 

As a concrete application of the results of §5, which also exhibits the main features of the 
general iV-dimensional case, we consider in this section the general form of the evolution 
equation (|5.20|) in two dimensions. Then, we show how evolution equations that are used in 
statistical physics as model kinetic equations can be derived in the present context, at the 
same time throwing some light into their foundations from a new perspective. 
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An example of this type is Kramers' equation which for simplicity is considered here in 
one spatial dimension x: 



d t p + yd xP + F(x)d y p - Pd y (yp) - h 2 y p = (7.1) 



This is a kinetic equation for the probability distribution p of a (unit mass) brownian particle 
of velocity y in an equilibrium bath of temperature h/2f3, f3 being its friction (drift) coefficient 
and F an external force field. 

Eq( |7.1| ) can be derived from microscopic dynamics, e.g. by expanding the so-called gen- 
eralized master equation in powers of the ratio of the masses of the bath particles to that of 



the brownian particle, for brevity taken here equal to 1 (see e.g. [Pq] , ch.IX.4, §5). We 
notice however, that ( |7.1| ) has been extrapolated far beyond the usual model of brownian 
motion, namely, a heavy particle interacting via hard collisions with a bath of much lighter 
particles (e.g. it has been used for self-gravitating systems, see the discussion in p| §5 and 
references therein). 

On the other hand, ( |7.1| ) is also obtained by using Langevin's equation as a dynamical 
model and assuming, that the motion of the particle is a Markov process, and that motion 
along x is deterministic so that the dependence of the transition probability on x is a 5- 



function (see e.g. 153] eq(242), 031 ch.II.2, [ 40 1 §VIII.7). This is an extra condition, the use 



of which is reflected in the derivation of ( [7.1[ ) in the context of StC by writing Langevin's 
equation as a stochastic differential equation 

dx = ydt dy = -((3y - F(x))dt + hdW t (7.2) 

where Wt is a standard Wiener process and no stochastic term appears in dx (see e.g. [ (47j 



theorem 9.3.1, |39fl , §5.3.6, ||j §11.3). This extra condition is often justified by saying that 



the external field varies slowly on a length scale in which the velocity is damped (e.g. [iOfl, 
p. 232). Nevertheless, there is no generally accepted view on how ( [7.1|) should be generalized if 
this extra condition is relaxed, although second order derivatives in x are expected to appear 
on the basis of a systematic analysis of the microscopic dynamics of realistic classical 3D- 
models of brownian particles weakly coupled to the bath (]PJ §2, JiSJ §4).[] Such 2nd order 



derivatives are not easily accounted for by approaches based on StC, since this presupposes 



the appearance of a stochastic term in the first of eqs( |7.2|) , not easily justified from a physical 
point of view. 

On the other hand, in §5 we have shown how such evolution equations (including 
( [7.ip are formally incorporated as hamiltonian equations in symplectic mechanics developed 
in the context of the noncommutative DC defined by (|5.21|) . Below we explore the physical 
interpretation of this formal result by employing the techniques and conceptual framework 
of the previous sections. 

In the notation of §5 we write for a (2+l)D-system 

{u") = (u,v,w) , (x*) = (t,x,y), ( ai ) = (a,a), P» = (p,q,r) , X = ^P"0 +u „ (7.3) 

A* 

6 The same conclusion follows from a similar analysis of quantum systems and by taking the classical limit 
of the resulting equations; e.g. for a harmonic oscillator in a bath of oscillators, such an analysis leads to a 
widely used kinetic equation of the Lindblad type (i.e. conserving density matrices), whose classical limit is 
a FP equation involving a drift and a diffusion term in the x axis ( |4^| § 3,4). 
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{At) 



Then, by ( pi] ), C ^15|) , ( fTT9l) , we have 

X(x) = (dx,X) = 
X( J /) = (dy J X) = 



111 

K X LI 

k' X' li' 



a(np + Xq + fir) = -j—R x 
n u 

a 2 

a(n'p + X'q + /i'r) = - — 
^22 



(7.4) 



(7.5) 
(7.6) 



= -t{k p + X q + lm r) 



,12 



r/<7 



— (kk'p + XX' q + Lifi'r) 



(7.7) 



V 
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2 i \/2 i 12 \ 

k p + X q + li r) 



where the limiting probabilities in the continuous limit 
( CT ) to be 



(7.8) 

are readily found via (17.41) , 



(p, q, r) 



Xfi' — X'fi n'fi — kfi' X'k — Xk' 



\A\ 



\A\ 



\A\ = det{A 



(7.9) 



Finally, the evolution equation ( |5.18| ) is 
^ 2 a 2 „ „ „ 1 



^ " ^T R ^f - TJT R v d *f ~ ^ (Vndlf + m2 d 2 xy f + V22 d 2 f) + o(a, a) = (7.10) 
van OH22 <Z 

Now suppose that (x, y) are the phase space coordinates of a particle, moving under friction 
linear in the velocity y, say — f3y, and an external field F(x). Evidently, the generator of 
newtonian motion is (y, — (f3y — F(x)) and therefore as in §5 (cf. ( |5.13j )) we may introduce 
the 1-forms 

ai = dx + ydt , a 2 = dy — (f3y — F(x))dt 



Consequently ( |5.14| ), giving the "infinitesimal" changes of (dx, dy) in time dt, becomes 
{ai , X) = 0, i = 1,2, hence 



(da; ,X) = by = bR x , (dy , X) = -b((3y - F(x)) = bR y 
and in the limit 

R x = R x = y, R y = Ry = -((3y - F(x)) 



(7.11) 



(7.12) 



If motion along x is deterministic in the continuous limit, that is, the concept of a trajectory 
in x is meaningful, then by the remark (3) in §5, this is equivalent to 



r/11 = (hence i)i 2 = 0), that is k 2 p + X q + fi 2 f = 



(7.13) 
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This corresponds to the absence of a stochastic term in ( [7.2ft ). Since p ,q,r are nonnegative, 
it is readily obtained from ( |7.13| ), (|7.9| ) that up to a permutation of the lattice coordinates 
(u,v,w) (hence of (p,q,r)), there are two possibilities: 

(1) q = f = 0, p = 1, k = k' = implying fj 22 = 0, A//' — \'(x ^ 0. Then, (|5.20| ), the evolution 
equation in the continuous limit, is 

dtf - ydj + (Py - F(x))d y f = (7.14) 

This describes motion with velocity damping and well defined trajectories in phase space. 
There is a 4-parameter gauge freedom for the transformation (|7.4j) . 

(2) 

p! k' , , 

<? = 0, P=— zj, r = — -, /c = /i = 0, T] 22 = -Kp >0 (7.15) 

pi — K K — pi 

with the evolution equation ( |5.20|) being 

dtf - ydj + (Py - F(x))d y f + \(h 22 klp!) d 2 f = (7.16) 

and there is again a 4-parameter gauge freedom subject to the constraint k'pl < 0. This is 
formally identical with the adjoint of ( |7. 1| ) . This formal similarity is more than a coincidence 
and gives another illustration of the consistency of the probabilistic framework introduced 
in §3 and elaborated in §5. A simple choice of the gauge is k' = —pi = 1, A' = 0. Then, in 
view of Q, fl5TT6|) , eqsfpJl), QQ3D imply 

a 2 _ 1 

rf 2 = —{p + r ) + o(a, a), -(dy • dy, X) = rf 2 + o(a, a) 

Therefore 

1 a 2 a 2 a 2 

-(dy • dy, X) = — (p + r) + o(a, a) = —(1 - q) + o(a, a) = — + o(a, a, q) -> f\ 22 = h 22 

i.e. a 2 /b is the average of the square of the velocity change from y to y + dy in a "small" 
time interval b divided by b. 

On the other hand, it is well known that in the usual interpretation of ( |7. 1| ) as a model 
of brownian motion, h = Py 2 , with y =thermal velocity, whereas, for short times t in which 
the velocity changes from y to y, ((y — yo) 2 )/t = Py 2 + o(t) = h + o(t) with (, ) denoting 
the average over the bath (see e.g. |24j] eq(161)). Thus, in the continuous limit h 22 = h and 
( [7.16Q coincides with the adjoint of (|7. 1|) . 

It should be noticed here, that this result follows by assuming that in the configuration 
space (i.e. in the x-axis), motion is along well-defined trajectories. It is a remarkable fact 
however, that if this assumption is relaxed, then ( |7. 10|) reduces to an equation describing 
diffusion in the well, a fact not easily accommodated in other approaches, as 

explained at the beginning of this section. 

Remarks: (a) In the present section and in §4, we have seen that up to a coordinate 
transformation, our approach allows for concrete results in specific cases, once the drift term 
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(dx l ,X) i.e. the average "infinitesimal" change of x % is given, that is, once a prescription is 
given for the choice of X in each particular case. In [13], X was specified by assuming it to 



be hamiltonian in the context of the "second order calculus" ( |5.21|) . This points out to the 
need of developing symplectic mechanics on the oriented hypercubic lattice, to be presented 
in another paper. 

(b) A similar, but simpler procedure yields in the ID case the following: We write (|5.1|) , 
>73l) as 



a 2 



t = —b(u + v ) , x = a(Ku + Xv) , > h 

b 

Then using ( [4.1|) , eq( |5.16|) becomes 



A k - 

p = — - , p = - , Xk < 

k — A k — A 

and the evolution equation in the continuous limit, eq( p.20| ), gives 

d t f-Rd x f + h(Xk)dlf = 

Notice that (dx, X) /b = a(np + Xq)/b = a 2 R/bh — > R so that the results of §4 are recovered 
if k = -A = 1 (cf. eqs(PD-(pHp). 

8 Comments on higher order equations 

The basic methodological "rules" of the present approach introduced in the previous sections 
are: 

• The extended phase-space has the structure of a (oriented hypercubic) lattice Ai. 

• The phase space coordinates (x M ) are obtained by a coordinate transformation from 
the lattice coordinates (m m ), involving the scaling parameters (a M ), so that the usual 
continuum description in terms of the (x^) is recovered in the limit (|5.3|) (cf. [32fl, 
p.244). 

In principle however, one can imagine other limits as well. In this section we first examine 
this possibility in a somewhat more general setting and then we reconsider the case of the 
oriented hypercubic lattice in the light of the results obtained. 
Let the DC on M. be defined by 

du^ • du v = C£ u du p , C7 G A (8.1) 

where the summation convention has been used (it is used throughout this section). This 
may be considered as a class of (algebraic) deformations of the ordinary DC. It can be shown 
that (@ §4.5) 

df = D p fdu", f = f(u) E A 
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+00 .. 



(8.2) 



r=2 



f*l,— ,l*r 
vi,...,v r _2 



We pass to the continuous limit by dividing {w^} into several groups, with the same scaling 
in each group. In each case one of the groups contains the "time" variable(s). It will be 
readily seen, that the results obtained are valid even if in each group, each coordinate is 
scaled differently, provided the order of magnitude is the same for all elements of the same 
group. Moreover, any finite limit of the scaling parameters is set equal to 1, to avoid useless 
cumbersome notation. 



(1) Division of into two groups: 

t a 

u a =: — , a — 0, . . . , n , 
b 

continuous limit: 



u 



x" 
a 



n + 1, 



N 



a, b — > , — 



From (O) and (O) we get 



(8.3) 



dt a • dt b 
dt a • dx i 
dx l • dx J 



bCf dt c + — Cf dx\ 
a 

aCl 1 dt b + bCf dx j , 
a 2 

— C^dt a + aC l k °dx k 
b 



•4) 



daf + \bC b a c d b d c f + a C« d b d t f + ~j C* didjf + o(a) 
^f + ~ Cf d a d b f + b C: j d a djf + \a (fOAf + o(a) 



D a f 
DJ 



If as in §5, quantities with a hat denote values in the limit (|8.3| ), then in that limit, 
reduce to (cf. eqs(Pfl, (^Tj) and @ eq.(4.13)) 



dt a • df = 0, 



dt a • dx l = 0, 



dx l • dx 3 = C%dt a 



■5) 



5.6) 



D a f = d a f+-C% didjf, Dif = dtf (8.7) 

so that the evolution equation (|3.34f ) reduces to a well-defined equation of the form (|5.20|) . 
(2) Division of into three groups: 



t a 

u a =: — , a = 0, . . . , m , u T 
b 



r = m + 1, 



x 



i = n + l,...,iV 
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continuous limit: 



a, 6, c 



Eqs (B, fj) imply 



dr • dt b 
dt a • dy r 
dt a • dx l 
dy r • dy s 
dy r • dx l 
dx l • dxi 



a 2 _ a 3 

°> T 

c o 



6 Cf dt c + - Cf dy r + — Cfdx 1 
c a 

DC 

c Ci r dt b + b C* r dy s + - CfW 
, a 

a Cf dt 6 + — Cfdy r + 6 Cf dx j 
c 

c 2 c 2 

— C 7 n s dt a + cC^dy 1 + - C[W 

— C™dt a + a C s "dy s + C qw 
y + — C?'dy r + a C*Mx fc 



•9) 



Dif = dif + o(a), 



Drf^drf+I-C^d^f + Oia 3 ) 

— c 



5.10) 



1 (7 r<7 



(8.11) 



It is now clear that in the limit flS.8|) , eqs (|8.10| ) and all but the last one of (|8.9| ) have a 
well-defined limit. The coefficient of C l J in the last of (|8.9|) and in (|8.11|) diverges as 1/a. 
Therefore, the limit (|8.8|) exists in this case only if we impose an extra condition on the 



structure constants which define the DC 

In this case (|8.9|) - (|8.11|) reduce in the continuous limit (|8.8|) : 

dt a • dt b = 0, dt a • dy r = 0, dt a • dx l = 0, dy r • dy s = 0, 

dy r • dx i = C r Jdt a , dx i • dx j = K'^dt a + C l r j dy r 



[8.12) 



(8.13) 



Dif = dif, D r f = d r f + ~C?d i d j f 

D a f = d a f + Cfdrdif + ^ d.dj f + ^C^djdJ 



(8.14) 



so that the evolution equation (|3.34|) reduces to an equation containing 3rd order derivatives. 
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It is easily seen that proceeding in this way, that is, by dividing into disjoint subsets, 
scaled by parameter a a, A = 1,2, ... ,L and then considering the limits 

a A ^0, &1^_> finite, CI? -> finite (8.15) 

a A p 

then unless L = 2, e:rfra conditions on Cj? should be imposed so that in that limit the 
commutation relations and df, f G A, are well-defined. 

The significance of this general result can be better appreciated if we implement the 
procedure described above, in the case of the oriented hypercubic lattice M., that is, when 
the DC is that of §5, ( |5.8|) . Instead of the limit (|5.3|) (or equivalently (|8.3|) ), we consider 

a„-0, WpL^hq* y-±oo, K^k (8-16) 

that is, essentially the limit ( |3.8|) (to simplify the presentation, no coordinates of order b 2 ^ 3 
are considered since they lead to the results of §5). However, it will become evident that the 
results obtained are valid for any limit of the type ( |8.15| ). 
By (|8.16| ) we have 



ai = bai, CLiCLjCLk -> hijk, b := b 1/3 (8.17) 

For arbitrary / G A, df is given by (|5.5|) and we consider its value in the limit (|8.16|) : 
Expanding ( |5.6|) and using that Y^n^i = Y^u^-^i = 0' once again we get (|8.10| i), namely 



Similarly, ( |5.7| ) gives 

+ l E ( E MJ) Mi) Mi) 5 o ) d^dif + o(6) (8.18) 
To simplify the notation, we put 

AXO := E Mte - 4(0, e := fe) G 

i 

and we consider the functions 

^(0 = ^T,(MZ)) 2 BZ, M0 = IJ2(MZ)) 3 b%, (8.19) 



related to the coefficients in (|8.18|) . We notice that by the Schwartz inequality 



1^(01 < nen iiaj < lien max{||Aj}=:M|ieii ^Miiei 

At 
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where A, 



(My 



and || || is the Euclidean norm in R . Therefore, 



\M()\ < (MO? \BS\ 



4(0 



-°0 



(8.20) 



where we used that Bq 



B% > 0, eq (^T7D. 



On the other hand, 

\2 



8.18|) is well defined in the limit ( |S.17|) , then 



if we require that 

^2 (^niOY Eft in ( |3.19| i) should be at least of order b (cf. ( |3.12| )) and therefore, its limiting 
value should be zero, i.e. 



R* 1 

n 



Therefore, $3(£) = for all £ G R , that is, the coefficients of 3rd derivatives of / in (8.18) 
vanish. 

We remark the following: 

(i) From the form of ( |3.18| ), QS.19 ), it is clear that if we require df to be well-defined in the 
limit (|8.17| ), then all coefficients of the derivatives of order higher than the second in the 
expansion of df, eq( |5.5| ), vanish. 

(ii) The same conclusion holds if we consider the limit ( |8.15|) , that is, instead of 



a { = b^ L a h L > 3. 

(iii) For the above result, it is essential that Bq > 0. This is a consequence of our interpre- 
tation of vector fields as transition probability distributions. In fact, the procedure followed 
above, is a modification of an argument used in the theory of Markov processes which ensures 
that higher than the second moments of such processes vanish ( ||39|| , §3.4). 

(iv) It can be seen from ( |3.18|) , that in the limit ( |3. 17j ), the resulting second order coefficients 



are no longer nonnegative-definite (see also (|8.14|) ). 

Therefore, summarizing our results in this section, we may say that a continuous de- 
scription can be obtained from any discrete lattice structure only if the discrete phase space 
coordinates are at least of order \fi, where t is the time. In this case we get evolution equa- 
tions at most of the 2nd order with nonnegative-definite 2nd order coefficient, which is the 
most general (linear autonomous) differential generator conserving probabilities (P0[ §§4,5). 



9 Discussion 

In this section we summarize the basic assumptions introduced in this paper, results and 
conceptual insights obtained and we comment on the limitations of the present approach 
and on directions in which it can be further elaborated. 

We started with a discrete picture of the extended phase space M.. This corresponds 
to looking at all fine structural details of the system and subsequently pass to its "coarser" 
description, which corresponds to the limiting procedures followed in §§4-8, as for instance 
done in statistical mechanics and lattice field theory. 

On the other hand, we saw in §3 that a discrete picture allows for a geometric visualiza- 
tion of the universal DC on M. as a di-graph, from which special DC on M. can be easily 
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constructed. In this way, their necessarily noncommutative character discussed in §3.1 is in- 
terpreted as due to the non- vanishing "size" of the differentials. In fact, these general ideas 
were illustrated in §2. The ID discrete model studied there, clearly suggested connections 
between motion defined on Ai endowed with a noncommutative DC, and a random walk 
on an oriented square lattice, one dimension of which is related to time (cf. ( |2.4| )). In §3, 
vector fields X were defined as elements of the space dual to 1-forms and turn to be 1st order 
difference operators on the algebra A of functions on Ai, when M. is an (N + l)-dimensional 
oriented hypercubic lattice. Here, a crucial fact is that only vector fields of a special form are 
the generators of automorphisms of A, in contrast to the usual (commutative) DC, where 
any vector field X (i.e. 1st order differential operator) generates (local) automorphisms of 
A. That a vector field generates automorphisms of A is equivalent to the fact that it induces 
a flow of well defined trajectories on M.. In the discrete case, these are paths which connect 
points of the lattice. This raises the question of what kind of motion is described by an ar- 
bitrary vector field Thus, we have been led to a new way of looking at vector fields, namely, 
both as generators of evolution of observables and as states describing transition probabilities 
for "infinitesimal" changes on the lattice. 

In this new perspective, probabilistic concepts are introduced in the dynamics at the 
"infinitesimal level", in contrast to (classical) statistical mechanics, where microscopic dy- 
namics and probability distributions are two a priori quite distinct concepts. Mathematically 
speaking, this double role of vector fields stems from the fact that in the present context, the 
concept of a vector field is much more general than the generator of an automorphism of 
A, as in the case of classical dynamics based on ordinary DC (cf. the discussion at the end 
of §3.2). In fact, this is a general feature of noncommutative DC on both discrete and "con- 
tinuous" manifolds, which allows for the description of evolution equations involving higher 
than the 1st derivatives and/or difference operators: By using the concept of a vector field, 
a noncommutative version of differential geometry and tensor calculus can be developed in 
close analogy with the ordinary (commutative) case. This is of potential value in many areas 
of physics (see e.g. pOfl , fl2| , |0§). We will come back to this point at the end of this section. 



An indication for the consistency of the interpretation of vector fields X as giving the 
transition probability distribution for "infinitesimal" motion along the lattice axes, is pro- 
vided by the proof that the associated correlation matrix P = (P^) for the differentials of 
the lattice coordinates, vanishes if and only if X generates automorphisms of A, or equiva- 
lently, if trajectories along the lattice are well defined. This means that there is no possible 
interference of motions along different axes at the same point (motion along a specific axis 
is either impossible, or certain). Moreover, in this way it became clear that, a function / 
on M. is an observable, in the sense that expectation values for its change represented by 
df, are given by (df,X), multiplication of observables corresponding to the •-product in 
O 1 , eq (|1.3| ). On the basis of this and assuming on physical grounds that time flows with 
certainty, i.e. for every vector field X there is always a change of some fixed element, t say, 
of A, it follows that dt is proportional to the (essentially) unique eigenvector of the above 
mentioned correlation matrix (P^) belonging to the zero eigenvalue. It turns out to be the 
unit p of the commutative algebra (^(A), •) of 1-forms, dt = —bp (§§3.3, 3.4). Using these 
results, the generator of dynamical evolution, X, is determined by iV 1-forms a 1 vanishing 
along X, eq( |5.14| ) and the time evolution of observables is given by —X(f)/b = 0, eq (|3.34|) , 
in direct analogy with classical dynamics. 
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In §§4, 5, using the probabilistic framework of §3, we considered appropriately scaled 
linear coordinate transformations on the lattice. We passed to the continuous limit in which 
M. becomes M> N+1 in the same way this is done in the theory of brownian motion. We showed 
that the evolution equation for observables is a 2nd order partial differential equation with 
nonnegative-definite leading coefficient. This is the limit of the correlation matrix (P M ^) 
of (changes along) the lattice axes, namely (fj^) in ( |5.19|) or (|5.25 ). Each of its elements 



measures the deviation of the evolution operator from an ordinary derivation (i.e. generator of 
deterministic motion) in the corresponding phase space axis. Thus, noncommutativity of the 
DC in the continuous limit (that is, in a "coarse-grained" picture of the system's evolution) 
is due to the fact that on a discrete ("fine") level, motion described by a vector field, is in 
general not along well defined phase-space trajectories. This is an idea also appearing in 
the theory of brownian motion and stochastic mechanics, though in a completely different 



mathematical and conceptual framework (see e.g. |48|], |]49 |, |42|| ). In fact, the explicit form 
of (fj l 3) in ( p.25| ) (cf. ( [5.23D as well) is consistent with the form of the diffusion matrix 
in the kinetic theory of open systems, the theory of markovian stochastic processes and 
StC. Moreover, the coefficients of 1st order derivatives in the evolution equation ( p.2(J| ), are 
just the first moments of the coordinate changes with respect to the transition probability 



distribution defined by the vector field X (cf. ( |5.11|) , (|5.15| )). Therefore, our evolution 



equation ( |5.20| ) is the formal adjoint of a generalized FP equation (backward equation), 
justifying our interpretation of / G A as observables and thus giving still another indication 
of the consistency of our approach. An interesting conclusion in this context is that the 
interpretation of vector fields as transition probability distributions on the lattice, implies in 
the continuous limit that evolution is forward in time. Thus, if we accept the correspondence 

discrete description < — > fine details of the system's microstructure 

i i 
continuous description < — > "coarse grained" (macroscopic) picture 

then, in the present conceptual framework, irreversible evolution is a characteristic of macro- 
scopic systems. 

The continuous limit considered in §5, at first sight may appear an arbitrary choice among 
many possible ones that would lead (via the expansion of ( |5.10| )) to evolution equations in 



general involving derivatives of any order. This issue was studied in §8 and we have shown 
that (i) the limiting procedure of §5 is the only way to pass to a continuous description by 
scaling the lattice coordinates, without imposing additional ad hoc conditions on the DC, 
that is, it is the only continuous limit independent of the general form of the DC one starts 
with; (ii) if the probabilistic interpretation of vector fields introduced in §3 is employed, then 
for all types of limiting procedures, higher than the second order differential equations are 
impossible as evolution equations. Thus, (i) and (ii) imply the unique character of the limit 
considered in §5. 

In §§4, 7 we considered simple applications of our approach: 

(a) The derivation of the ID diffusion and Smoluchowski equations for a constant external 
field, and the ID FP equation in velocity space giving the Ornstein-Uhlenbeck process, thus 
illustrating the possibility to incorporate random walk models in the present context. 

(b) The derivation of Kramers' equation in one spatial dimension, by assuming (i) (new- 
tonian) motion under friction linear in the velocity and an external field; (ii) trajectories 
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exist in configuration space, that is, motion is deterministic there. In fact, by using the 
explicit form of the correlation matrix ffi under the above assumptions, we have shown that 
the diffusion coefficient in our evolution equation coincides with the diffusion coefficient of 
Kramers' equation computed in the context of the theory of brownian motion, or of kinetic 
theory. 

On the other hand, our result makes explicit the use of assumption (ii) above, a fact often 
hidden in the derivations of Kramers' equation from physically plausible stochastic models. 
Although this is also made explicit when Kramers' equation is derived from Langevin's 
equation seen as a stochastic differential equation, in our opinion, the present approach 
has some advantage: In the context of StC, it is not easy to relax this assumption (see 
the discussion following (|7.2| )). In the present context however, this is possible. Actually, 
the general form of the corresponding corrections is evident, namely, a diffusion term in 
the configuration space, much in accordance with what seems plausible on the basis of the 
(classical and quantum) kinetic theory of open systems. Finally, it is clear that these results 
are valid in the 3D case. 

The approach elaborated in this paper, clarifies several conceptual issues concerning the 
relevance of noncommutative DC to (the derivation of) kinetic equations and to StC. In 
addition, by giving a definite prescription for the general form of irreversible evolution equa- 
tions in terms of vector fields X in the discrete framework, it provides an adequate discrete 
formalism for deriving such equations in an appropriate continuous limit. However, it should 
be supplemented by a "dynamics" , that is, a general procedure for choosing the generator 
X in particular cases. In sections §5 and 7, we have used Newton's equations written as 
1-form relations (see ( f5.13j ) and the derivation of ( |5.15| ), ( |7.12| )). As briefly discussed at the 
end of §7, a more systematic approach is to develop symplectic mechanics in the present 
context and require X to be hamiltonian, in analogy with classical mechanics. In fact, from 
a mathematical point of view, this approach has been followed successfully in the continuous 
regime for the second order DC defined by ( |5.21|) , giving the promising result that hamilto- 
nian evolution equations have the form of generalized FP equations that appear in kinetic 
theory ([0). In another paper, the present approach will be elaborated in this direction. 

It is also possible (and physically desirable) to consider the extension of the present 
formalism when the transformation from the lattice coordinates (u^) to the scaled coordinates 
(t,x l ), eq(p.l[), is not linear. Many of the present results are expected to be still valid. 

Finally, the general formalism and the conceptual framework in this paper, may be ex- 
tended from the case of the oriented hypercubic lattice, to more general structures induced 
by appropriate DC on a discrete manifold. More precisely, it is clear that throughout this 
work, differentials are not necessarily infinitesimal quantities in the usual geometric sense of 
ordinary DC. Rather than that, if seen in the suggestive representation of a DC as a di-graph, 
they express interrelations between points of the discrete manifold, so that all points related 
to a given one, should be considered as being neighboring to it (see (|3.5|), its interpretation 
in §3.1 and the differentials for the hypercubic lattice, eq( [3.21|) ). This idea, that a DC is 
based on a concept of "relational" , rather than geometric type of "proximity" may be further 
elaborated. In this way, it may become possible to develop a general mathematical and con- 
ceptual framework for describing physical systems based on the interrelations (interactions) 
among its different parts, rather than on their relative geometric position. This is virtually 
relevant in situations where collective effects are significant, or even dominant (e.g. systems 
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with long-range interactions), for which traditional approaches often do not work beyond 
the lowest approximation. We will come back to this point in another work. 



Appendix A 

1. Here we prove the following proposition referred to in §3.2. 

Proposition A 1-1 and onto mapping $ : M. i— ► M. induces an automorphism <ft : A i— > A 
that maps the basis {e, , i E M.} of A onto itself. The converse is also true. 

Proof: : We define 

<j>:A^A: = /(*(«)) (A.l) 

Clearly is an endomorphism and 

0( e i)(j) = e i($(i)) = = 8*-H%)j = e*-i(i)0") 

hence 

0(e;) = e$-i ( j) (A. 2) 

which shows that is 1-1 and onto. 

"<£=": Conversely, let be an automorphism of A From (j3~2"li) we get 0(ej)(0(e.,) — <%) = 
0. If 0(ej) = ^2 k fik^k, the above equation implies that = or 1 , V fc. By (|3l2"|ii) 
^2ikfik e k = J2k e k hence Y^ifik — 1 where we have used that 0(1) = 1 and {ej} is linearly 
independent. Consequently, in each column k of fa there exists exactly one element, on the 
ifc-th row, say, which is nonzero and therefore, necessarily equal to 1. Now suppose that for 
k ^ I, the corresponding nonzero elements are on the same i-th row. Then, these columns 
are identical which is impossible, since is 1-1 and onto. Therefore, in each row i of fa, 
exactly one k is nonzero and equal to 1. Hence, 0(e«) = e& for some k depending on i. Now, 
0~ 1 ( e «) = e Hi) defines $ : M. i — > M. which is 1-1 and onto and satisfies ( |A.1| ). QED 

2. By considering the case of classical dynamics, we explain below the choice of the sign in 
( |3.33| ) and comment on the form of (|5.14|) . 

Let the equations of motion of an iV-dimensional dynamical system be 

The time evolution of observables / and states a, is given by semigroups of Koopman and 
Perron-Frobenius operators respectively, whose infinitesimal generators lead to 

d f J*, ^/ da d , { v 

with the summation convention used here. The operators on the r.h.s. are formal adjoints 
to each other (cf. e.g. [E3] ch.7). 
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Then, the prescription for extended dynamics is 



. . dx l . dt 

equations of motion: — — = H (x) , — = 1 

as as 

(d d s \ 

evolution of states: — + -— R l (x)- )a =: X s (a) = (A. 3) 

\ot ox 1 J 

( d d \ 

evolution of observables: — — + KUx)—^ if =: Xo(f) = 

V at ox 1 J 



Therefore, the system of characteristic equations of the evolution equation for observables, 
is given by the vector field with components (X (t), X (x 1 )) = (— i.e. evolution is 
obtained by integrating the equations of motion "backwards in time".Q This explains the 
(— ) sign in ( |3.33| ) and the form of ( |5.13| ). In fact, to the system of characteristic equations 
of the observables' evolution equation (|A.3| iii) are associated the 1-forms 

a* = dx* + iTdt z = l,2,...,iV (A.4) 

vanishing along the integral curves of the vector field Xq- This is equivalent to 

(a\X o ) = i = l,2,...,N (A.5) 

( , ) denoting the contraction of 1-forms with vector fields in the ordinary DC. In the con- 
ceptual framework of the present paper, this is given by ( |5.13| ), Q5.14Q . It is a completely 
general formulation that includes as special cases (in an obvious notation) 
Newton's equations (N = 2M) 

(3 r = dq r + v r dt , 7 r = dv r + fdt 

or Hamilton's equations 

an f)fJ 

^ = dq r + —dt, lr = d Pr -—dt 

OPr Oq r 

To connect this with the usual formulation of symplectic mechanics, we notice that, by 
introducing the symplectic matrix 

(^) = W=(_°j o 

(/ being the M x M identity matrix), the 1-forms above can be rewritten in a form identical 
to ( |5.13| ), namely 

air air 

a' = dx { + V J ij — dt ^ on = V - ^— dt 

4^ dxi ^ 3 dx l 



7 cf. the theory of stochastic processes and stochastic differential equations (e.g. |f7j], §2.6), where, for 
observables (i.e. bounded functions of a stochastic process): Generator of the evolution of observables = 
(backward generator of evolution of the probability distribution) = — (adjoint of forward generator of evolution 
of the probability distribution). 
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(i — 1,2, ... , 2M) so that the «j can be given in terms of the symplectic form 



uj := — Jijdx 1 A dx^ — dH A dt, a { - 

i,j 

where i is the interior product operator for ordinary differential forms. 
3. Some of the results in §3 can be seen in the light of the discussion in the previous 
subsection: By fl3.24|) , X{u^) = P M is the defining relation for the components of an arbitrary 
vector field X. The approach in this paper leads to concrete results in particular cases, once 
X is chosen appropriately in each case, say (du 11 , X) = : P M = as it was done for instance 
in §4. This fact can be recast in a more suggestive form, namely 

(5", X) = 0, where a M := du^ - Rfp (A.6) 

that is, the (N + 1) 1-forms <S M vanish along X. Notice that, owing to ^ M -P M = 1, we must 
impose the condition ^ M -R M = 1, hence J^c^ = 0, that is, the <S M are linearly dependent. 
Moreover, by ( [3.32D , the correlation matrix of the lattice coordinates is P^ v = (a' 1 • a v ,X) 
in close analogy to (|5.24|) . In fact, ( |A.6| ii) is just ( |5.13 ) written in the lattice coordinates 



To see this, we transform ( |A.6| ) to {t,x l ) coordinates. Since P M = R^, the become 
ciiA 1 ^ = Mj,du" - a i Ai ^P = dxi + ^(AR) l dt 



By setting (ai/b)(AR) 1 = R\ the second equation becomes (|5.13| ), hence ( |A.6| i) is transformed 
to (Ell. 



APPENDIX B 

As mentioned in §5 remark (5), in specific cases, instead of expanding the evolution equation 
( |5.10| ) as in (|5.12|) , it is possible to write it explicitly in terms of appropriate difference 
operators which in the continuous limit gives the corresponding form of ( |5.2U| ). As an 



illustration, we present here an iV-dimensional example which generalizes that of §2. We 
employ the definitions and notations of §§3, 5. 

The transformation ( [5.1|) is chosen as follows (i = 1,2,..., N) 



t = -&^V, x i = a i (u° + --- + u i - 1 -u i + u i+1 + --- + u N ). (B.l) 

so that its inverse ( |5.4| )and the commutation relations ( |5.8| ) become 

n n - 2 t 1 A x i , 1 ft x l \ . . 
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a 2 

dt • dt — —b dt, dt • dx % = —b dx\ dx l • dx l = —-^dt, 

b 

dx l • dx j = — C hp- dt + cij dx l + at dx j 



(B-3) 



b 

In the limit ( |5.3| ), eq( |5.21|) takes the form 



dt • dt = 0, dt • dx* = 0, dx l • dx j = -h^ dt, (B.4) 
For the differential of a function of 

d/ = • • • , u»-\ + 1, ...,u n )- f(u°, u n )) du». 

As a consequence eq (|5.5|) takes the form 

N 

- f(t - b, x + a, x l - Of, x + a) + (n - 2)f(t + b,x + a) + 2f(t, x) dt 



+ ^(dj) (t-b,x + a, x\ x + a) dx\ (B.5) 

with 

Ch/n <"....) : ^-[/( /'' • a,....) - j{ r'- a, ....)[ 

After a lengthy calculation we obtain 

N 2 

df = [(d^ t f)(t,x)-J2^-A^m-b,x + a,x\x + a) 



2b 

8=1 



(d +i d + jf)(t — 6, x + a, x 4 , x, x J , x + a)] dt 

l<i<j<N 

^2(Bif)(t - b, x + a,x\x + a) dx\ (B.6) 



where 



(d- t f)(t,x) :=±[f(t,x)-f(t-b,x)} 



(A,/)( r'....) :=![/(..., x l + a,,... ) + /(..., x^- a,,... )-2/(...,x i ,...)]. 



a? 1 



(cf. (l2l)-pT0|)). Evidently, in the limit (0), eq([E[6|) gives 

v 

d/=[(9J)(t,x)-^^ (#/)(*,*)+ 

i=l 

£ MW)(*> + Z) (B.7) 
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By ( |3.24| ), a vector field in this calculus has the form 



N N 



(Xf)(u) = J^P^+Mu) = ^"(/(u.u" + !,«)- /(«)) 




or, by dKT]) 



N „2 



(Xf)(t,x) 



X% x) [(3_ t /) (t, x) - ]T -i ( A J) (t-b,x + a,x\x + a) 



+ ^ (d +i d +j f)(t - b, x + a, x\x,x j , x + a)] 



l<i<j<N 



+ S ^2 / X t (t, x) (dif)(t - b, x + a,x l ,x + a), 



(B.8) 



X 1 = -b(P° + --- + P n ), X i = ai (P° + ■ ■■ + P^ 1 - P l + P i+1 + ■■■ + P n ) (B.9) 



and has + 1 solutions with P^ = 1, P v = 0, if v ^ fi, for fi = 0, 1, . . . , N. Using ( [EOD we 
readily get that in (t, x l ) coordinates, the evolution generator for observables, eq( |3.34j ), has 
components (l,±aj/6), thus obtaining the A^-dimensional generalization of the "random- 
walk model" of §2 (see eqflDjD). 
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